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Transverse momentum spectra, d 2 a / (dr/dp ^), of charged hadron production in pp-collisions are 
considered in terms of a recently introduced two component model. The shapes of the particle 
distributions vary as a function of c.m.s. energy in the collision and the measured pseudorapidity 
interval. In order to extract predictions on the double-differential cross-sections d 2 cr/(dr/dp^) of 
hadron production for future LHC-measurements the different sets of available experimental data 
have been used in this study. 


I. INTRODUCTION 

Recently a qualitative model considering two sources 
of hadroproduction has been introduced [I] . It was sug¬ 
gested to parametrize charged particle spectra by a sum 
of an exponential (Boltzmann-like) and a power-law pp 
distributions: 

V 7 2 ( ^Tkin/^e) T ~2 ) (1) 

d T} d Pj, (1 + t^n) N 

where Epkin = \Jvt + M 2 — M with M equal to the pro¬ 
duced hadron mass. A e , A,T e ,T, N are the free param¬ 
eters to be determined by fit to the data. The detailed 
arguments for this particular choice are given in [I] . The 
exponential term in this model is associated with ther- 
malized production of hadrons by valence quarks and a 
quark-gluon cloud coupled to them. While the power- 
law term is related to the mini-jet fragmentation of the 
virtual partons (pomeron in pQCD) exchanged between 
two colliding partonic systems. 

A typical charged particle spectrum as a function of 
transverse momentum, fitted with this function (|T]) is 
shown in figure |T] As one can see, the exponential term 
dominates the particle spectrum at low pp values. 

This model has already been shown to predict the 
mean < pp > values as a function of multiplicity in a 
collision [3j and pseudorapidity distributions of charged 
particles [1]. However, the major interest of many stud¬ 
ies in QCD is the transverse momentum spectrum itself. 
Therefore, in this article it is discussed how its shape 


* alexandr.bylinkin@cern.ch 

• nadezda.chernyavskaya@cern.ch 
1 rostov@itep.ru 



FIG. 1. Charge particle differential cross section 

(l/2irpT)d 2 a/dpdpT 0b] fitted to the function ([!]): the red 
(dashed) line shows the exponential term and the green (solid) 
one - the power law. 


varies in different experiments under various conditions. 
In [T] it was shown that the parameters of the fit (|T|) 
show a strong dependence on the collision energy. Un¬ 
fortunately, due to the fact that different collaborations 
measure charged particle production in their own phase 
space and under various experimental setup, the depen¬ 
dences observed in [T] were smeared and did not allow to 
make strong predictions for further measurements. Thus, 
an approach to correct the measurements in order to al¬ 
low an accurate combination of different experimental 
data is proposed here. 
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II. PARAMETER VARIATIONS 

In [4j it was shown that two sources of hadroproduc- 
tion described above, contribute to different pseudora¬ 
pidity regions: while the power-law term of 0 prevails 
in the most central region (77 ~ 0 ), the exponential term 
dominates at high values of 77 . Since each collaboration 
presents measurements on transverse momentum spec¬ 
tra in various pseudorapidity intervals, these variations 
might explain the smearing of the dependences in [I]. 
The idea to study parameter variations as a function of 
both collision energy and pseudorapidity region has al¬ 
ready been successfully tested in [5]. 

To further study the variations of the spectra shape as 
a function of pseudorapidity we use the data published by 
the UAl experiment [2] which are presented by charged 
particle spectra in five pseudorapidity bins, covering the 
total rapidity interval |? 7 | < 3.0. Figure [ 2 ] shows how the 
parameter N varies with pseudorapidity together with a 
power-law fit of this variation. Note, that the parameter 
shows a growth with pseudorapidity [ 3 ] , that is explained 
by higher thermalization of the spectra, as found in [4] . 



FIG. 3. Variation of the N parameter of |l| obtained from 
the fits to the experimental data na m no] (full points) as a 
function of c.m.s. energy y/s in a collision. Solid line shows a 
fit of this variation. In addition, open points show the param¬ 
eters for the data measured in another pseudorapidity inter¬ 
val 11 Tl 112) and pointed line shows the predictions calculated 
according to 



FIG. 2. Variation of the N parameter of (JTJ) obtained from 
the fits to the experimental data [2] as a function of pseudo¬ 
rapidity. Line shows a power-law fit of this variation. 

This variation can be parametrized in the following 
way: 

N = N°- (1 + 0.06- M 1 ' 52 ) (2) 

where N° denotes the parameter value at 77 ~ 0 and 
77 might be taken as a mean < 77 > of the measured 
pseudorapidity interval. 

Since the variations of the parameter as function of 
pseudorapidity have been found, it is desirable to ex¬ 
clude its influence when studying the dependences of N 
on the c.m.s. energy in a collision. This is possible, if 
one combine only those data that have been measured in 
more or less the same pseudorapidity intervals. Hence, 
we suggest to look first of all at ISR )B], PHENIX p] and 
ALICE [TO] data that were measured in the most central 
(|? 7 | < 0 . 8 ) pseudorapidity region. 


Figure [3] shows the parameter N variation as a func¬ 
tion of c.m.s. energy in a collision. One can notice that 
it can be characterized by the falling V-value. It is re¬ 
lated to the fact that the probability to produce a high-p^ 
mini-jet should grow with y/s. Notably, this behaviour 
correlates with the fact that N decreases when the rapid¬ 
ity interval between the secondary hadron and the initial 
proton increases. 

One can extract the following regularity from the fit 
shown in figure 


N = 2 + 5.25 • s 


- 0.093 


( 3 ) 


Remarkably, in the s —» 00 limit N —> 2, which corre¬ 
sponds to d 2 a/dpT oc 1 /l 4 i n the proposed parametriza- 
tion Such behaviour can be expected just from the 
dimensional counting in pQCD. However, in a real colli¬ 
sion one should take the initial conditions and the kine¬ 
matic restrictions into account, resulting in a higher value 
of N for lower y/s. Obviously, such initial conditions 
should become negligible in the s —> 00 limit, that is 
confirmed by the observed behavior (13]). 

In addition, figure [3] shows UAl pLl] and CMS [L2] 
data measured under different experimental conditions. 
In these measurements the pseudorapidity interval was 
much wider (I 77 I < 2.5) than in [ 6 ) 0 dD]. Therefore, 
one can compare the parameter values obtained from the 
fit of these data (open points in figure [3j) to the values 
extrapolated from Q with N° calculated according to 
([3]) and |? 7 | = 1.25 (pointed lines) and see a rather good 
agreement. Notably, this behaviour correlates with the 


1 In the and below y/s is given in GeV. 
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fact that N decreases when the rapidity interval between 
the secondary hadron and the same side beam proton 
increases. 

Let us now check this correlation explicitly and calcu¬ 
late the rapidity interval in the moving proton rest frame 
according to a simple formula: 

if = \v\ ~ log(^s/2m p ), (4) 

where m p is the mass of the incoming proton. The results 
of such a procedure are shown in the figure [4] Surpris¬ 
ingly, all the points came to a single line in this inter¬ 
pretation. To understand the origin of this universality 
one might use MC-generators: hard processes at large 
Pt are known to be described by MC generators pretty 
well, thus it is expected to get the value of N-parameter 
from the fits of the MC-generated spectra rather close 
to the real data, but with a higher accuracy and in a 
wider collision energy range. To check this universality, 
we have produced the Monte Carlo samples for proton- 
proton collisions at different energies for inelastic(INEL) 
events with the PYTHIA 8.2 generator [7], Indeed, the 
values of the parameter N extracted from the fits to the 
MC-generated spectra are nicely placed at the same line. 
Thus, a universal parameter describing the shape of the 
transverse momentum spectra in pp-collisions has been 
found. 


N = 5.04 + 0 . 2777 ' (5) 



FIG. 4. The dependence of the N parameter on the pseudo¬ 
rapidity of the secondary hadron in the moving proton rest 
frame. The data points from different experiments [51EIIH1E2- 
H2] are shown together with the generated data(MC). 

The further check of the observed phenomena can be 
made if one plot A as a function of the logarithm of 
the maximal kinematically allowed transverse momen¬ 
tum PT max of the secondary hadron at the specific pseu¬ 
dorapidity p. 


This dependence is shown in the figure [5] One can see, 
that larger N corresponds to the smaller PT ma . x - That 
should correspond to the x —> 1 limit of PDFs, where 
the fall-off of the perturbative cross section is modified 
by the fall-off of the PDFs. 



FIG. 5. The dependence of the N parameter on the maximal 
kinematically allowed transverse momentum of the secondary 
hadron. The data points from different experiments 
ll0lI12j are shown together with the generated data(MC). 


Remarkably, similarly to N, T and T e also show depen¬ 
dences as a function of both the collisions energy y/s and 
the measured pseudorapidity interval p. The variations 
of the T and T e parameters were studied in [5], In [5] 
the possible theoretical explanation of the thermalized 
particle production was presented and the following pro¬ 
portionalities were established: 

T = 409 • (v's) 0 ' 06 • exp(0.06|77|) MeV (7) 

T e = 98 • (Vi) 0 ' 06 • exp(0.06M) MeV (8) 


The parametrizations for T and T e differes only by a 
constant factor. However, both T and T e parameters 
reflect the thermalization which is stronger at higher en¬ 
ergies and when closer to the valence quarks. Therefore, 
the 081 parametrizations which are functions of centre 
of mass energy and rapidity interval can be rewritten 
in the form with only one universal parameter. This 
universal parameter is the rapidity distance p" from the 
farther incoming proton. 


p" = \p\ + log(Vs/2m p ) 

Using (7]8) we get the universal dependency 
T = 409 • exp(0.06j/') • (2 m p f m MeV, 

T = 98 • exp(0.06r/') • (2 m p ) 0 06 MeV. 
The dependences are shown on figure [6j 


(9) 

( 10 ) 

( 11 ) 


\fs 

PTmax = — sin[2 ■ tan -1 (exp(-\p\))] 


( 6 ) 


In the ( I0|ll|22 I m p is given in units of 1 GeV/c 2 
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FIG. 6. The dependence of the T and T e parameter on the 
pseudorapidity of the secondary hadron in the moving oppo¬ 
site side proton rest frame. The data points from different 

experiments EH5HSE3E2I are shown. 


III. PREDICTION FOR FURTHER 
MEASUREMENTS 

Though MC was shown to provide a nice description 
of the high-pt part of the spectra, the nature of the soft 
particle production still remains ambiguous and varies 
for different MC-generators. Therefore, the next step in 
understanding the underlying dynamics of high energy 
hadronic processes was done in a recent analysis [5]. 

In 0] it was shown, that the introduced approach is 
able to give predictions on the pseudorapidity distribu¬ 
tions in high energy collisions for non-single diffractive 
events (NSD). Using the parameterizations from 0] in 
addition with one can provide a formula that 

describes the shapes of charged particles spectra, being 
function only of a centre of mass energy and a measured 
pseudorapidity region. Let us now summarize all the 
equations to obtain the final result^ 


dcr \ a ( r 


dv) r ^~ Apower ' exp \ 2 aL 


power 


( 12 ) 


da 

di] 


— ) = A exp ■ exp 


exp 


A 

s*-exn 


■exp ■ exp 


(V ~ Vexp) 2 

2(J exp 

(V + Vexp) 2 


2 V exp 


(13) 


&power = 0.217 + 0.235 • In •> /s, 


(14) 


Vexp = 0.692 + 0.293 • In y/s, 

(15) 

<?exp = 0.896 + 0.136 • ln\/s, 

(16) 

A _ n i q . (.0.175 

power c/.xu o 5 

(17) 

A exp = 0.76 • s 0106 , 

(18) 


Now, with the knowledge of the variations of T, N 
and T e parameters and the exponential and power-law 
contributions one can calculate the normalisation pa¬ 
rameters A and A e in ([I]) in the following way. 


da 
dr] 


/ da \ / da\ 

V d V ) exp \ dr l)pow 



d 2 a 
drjdP'j, 


dprp 


(19) 




(1 + 


A 


Pt \n 

T 2 -N > 


dp?p = 


ANT 2 
N - 1 


( 20 ) 




A e exp> {-E Tkin /T e )dp 2 T = 2 A e T e (m+T e ) 


( 21 ) 

Thus, we get the set of equations allowing us to make 
a prediction for a double differential cross section, using 
the formula ([!]): 


\v\ = 


+ M min )/2, 


N = 5.04 + 0.2777' 

T = 409 • exp(0.06??") • (2 to p ) 0 06 , 
T = 98 • exp(0.0677") • (2 m p ) 

= 


0.06 


( 22 ) 


_ (JV-l) 

rr — NT 2 


2T e (m+T e ) (da/dr]) 
( da/dr]) pow , 


exp 5 


where r/ and rf’can be calculated using Q and 
correspondingly. 

Now, one can calculate double differential cross sec¬ 
tions d 2 a/ ( drjdpj,) of charged particle production in high 
energy collisions at different energies for NSD events. 
These predictions are shown in figure [7] for \r]\ < 0.8 and 
|? 7 | < 2.4 pseudorapidity intervals together with the ex¬ 
perimental data measured by ALICE JfJj and CMS f!2] , 
A good agreement of the prediction with the data can be 
observed. Thus, these results (22) give us a powerful tool 


for predicting the spectral shapes in NSD events. 


IV. CONCLUSION 


3 In the ( |l2| | and below A and A f: is measured in mb/GW 2 , A pow 
and A exp in mb 


In conclusion, transverse momentum spectra in pp- 
collisions have been considered using a two component 
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FIG. 7. Predictions the yield of charged particles 
(l/2npT)d 2 N/(dr/dpr) in high energy collisions in NSD events 
together with data points from ALICE [5] and CMS m ex¬ 
periments. 


model. Variations of the parameters obtained from the 
fit have been studied as a function of pseudorapidity ?] 
and c.m.s. energy yfs in the collision. A universal pa¬ 
rameter describing a shape of the spectra in pp-collisions 
was found to be a preudorapidity of a secondary hadron 
in moving proton rest frame. Finally, the observed 
dependences, together with previous investigations al¬ 
lowed to make predictions on double differential spectra 
d 1 2 3 4 5 6 a / (dp^dri) at higher energies, successfully tested on 
the available experimental data. 
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